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Abstract. In this work we present first results of our current project aimed at combining 
the 3D hydrodynamical stellar atmosphere approach with non-LTE (NLTE) spectral line 
synthesis for a number of key chemical species. We carried out a full 3D-NLTE spectrum 
synthesis of the oxygen IR 777 nm triplet, using a modified and improved version of our 
NLTE3D package to calculate departure coefficients for the atomic levels of oxygen in a 
C0 5 B0LD 3D hydrodynamical solar model atmosphere. Spectral line synthesis was subse- 
quently performed with the Linfor3D code. In agreement with previous studies, we find 
that the lines of the oxygen triplet produce deeper cores under NLTE conditions, due to the 
^ diminished line source function in the line forming region. This means that the solar oxygen 

IR 777 nm lines should be stronger in NLTE, leading to negative 3D NLTE-LTE abundance 
corrections. Qualitatively this result would support previous claims for a relatively low solar 
oxygen abundance. Finally, we outline several further steps that need to be taken in order to 
improve the physical realism and numerical accuracy of our current 3D-NLTE calculations. 
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1. Introduction by the realism of the ingredients used in the 

abundance analysis: atomic data that describe 
Photospheric abundances of chemical elements the spectral lines themselves, stellar model at- 
derived from stellar spectra are important mospheres, and the allowance of possible de- 
means for testing and constraining models of partures from local thermodynamic equilib- 
the formation and evolution of the Galaxy and rium (LTE). 

its various stellar populations. The reliability During the recent years, reliable spec- 
of the derived chemical abundances is limited, tral line data has been measured and/or corn- 
apart fix)m the quality of the observational data, puted for a large number of astrophysically 
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relevant spectral lines, and departures from 
LTE are currently accounted for with an in- 
creasingly more realistic treatment of colli- 
sional processes. However, the stellar model 
atmospheres that are commonly used for the 
abundance analysis are based on a number 
of critical simplifications: they are typically 
constructed in one-dimensional (ID) geome- 
try - plane-parallel or spherically-symmetric - 
and are subject to hydrostatic and radiative- 
convective equilibrium (see, e.g., Castelli 
& Kurucz 120041; iBrott & Hauschildtl 120051: 



Gustafsso n et al. 2008 ). Even though the input 
microphysics (opacities, equation of state) of 
these models is sufficiently realistic, the con- 
vective energy transport, one of the main pro- 
cesses shaping the physical structure of the 
photosphere, is treated in the approximation 
of the mixing-length theory (Bohm-Vitense 
1958) or its derivatives (ICanuto & Mazzitellil 
1991). 



In this respect, three-dimensional (3D) hy- 
drodynamical stellar model atmospheres are 
much more realistic, as they are able to account 
for convection based on first princi ples with- 
out the need of free parameters (see iNor dlund 



1982; Stein & Nordlund 1998; Asplund et al. 
2000; Frevtag et al.l 12.002c IWedemever et al.l 



2004). Moreover, this type of model atmo- 



spheres naturally allows for the emergence of 
a surface granulation pattern, horizontal inho- 
mogeneities and wave activity, all being unique 
properties of 3D hydrodynamical model atmo- 
spheres. It has been found that the interplay 
of these physical processes significantly alters 
spectral line formation, which may lead to sub- 
stantial differences in the abundances derived 
with 3 D and ID model atmospheres, respec- 
tively dCollet et al.l 120071: IDobrovolskas et all 



2010; Ivanauskas et al. 2010; Kucinskas et al. 
20131 : IDobrovolskas et al.l 12.01 3h . These find- 



ings conspicuously indicate that the 3D pho- 
tospheric structure and dynamics has a signifi- 
cant impact on the spectral line formation, and 
should be properly taken into account in the 
chemical abundance analysis. 

Besides the application of realistic stel- 
lar model atmospheres, an adequate treatment 
of non-LTE (NLTE) processes can be also 
very important for deriving reliable chemi- 



cal abundances. In the optically thin line- 
forming regions, the absorbing particles ex- 
perience a radiation field that is of non-local 
origin since it forms deeper in the photo- 
sphere. Consequently, radiation field in the 
line-forming region may exhibit significant 
deviations from the local Planck function. 
Radiation of non-local origin can distort the 
collisional ionization balance given by the 
Saha equation, and drive the population num- 
bers of the upper and lower atomic levels of 
the given transition away from the Boltzmann 
distribution that is valid in LTE (for a thorough 
discu ssion of these effects see, e.g. lBruls et al.1 
1992). Departures from LTE modify both the 
strength and shape of the spectral lines, and 
therefore can significantly alter spectroscopi- 
cally derived chemical abun dances, in particu- 
lar at low metallicities (e.g jThevenin & Idiartl 
Il999h . 

Given the significance of NLTE effects and 
the magnitude of the LTE 3D- ID abundance 
corrections (especially in the metal-poor stars), 
it is obvious that the two factors should be si- 
multaneously taken into account in order to 
derive reliable abundance estimates. However, 
the joint 3D-NLTE approach has only been ap- 
plied in a few selected cases so far (Asplund 
et al. |2004l:IShchukina et al.ll2005l:ICavrel et al.1 



2007; Pereira et al. 2009: ISbordone etal.1201 
Steffe n et al.ll2012h . therefore a systematic ap- 
plication of this admittedly very demanding 
methodology to, e.g., the investigation of stel- 
lar populations, is yet to come. 

Driven by the need to combine modern 3D 
stellar model atmospheres and non-LTE spec- 
tral line synthesis, we have recently started a 
project to make this methodology available for 
a wider range of chemical elements to be stud- 
ied with 3D hydrodynamical C0 5 B0LD stellar 
model atmospheres. Oxygen is a particularly 
interesting element to investigate using the full 
3D NLTE approach: it is the most abundant 
chemical element besides hydrogen and he- 
lium, and its photospheric abundance is widely 
used to trace the formation and chemical evolu- 
tion of various Galactic populations. It is well 
known that the O i IR triplet (A = 777 nm) ex- 
periences significant departures from the LTE 
(lKiselmanlfl993l : lFabbian et alJl2009L and ref- 



284 



Prakapavicius et al.: Oxygen spectral line synthesis in 3D-NLTE 



erences therein). Due to its high ionization 
potential, oxygen does not experience signif- 
icant overionization in late-type stars, so that 
NLTE effects are mainly limited to deviations 
from the thermal excitation equilibrium. These 
physical aspects limit the range of possible 
NLTE effects in the case of oxygen and there- 
fore make it a good test case for the full 3D- 
NLTE approach with C0 5 BOLD model atmo- 
spheres. 

2. Methodology 

2.1. Model atmospheres 

The 3D hydrodynamical solar model atmo- 
sphere used in this work was compute d with 
the CO 5 BOLD code dFrevtag et al.ll2012l) . The 
C0 5 BOLD code solves time-dependent equa- 
tions of compressible hydrodynamics and ra- 
diation transfer on a 3D Cartesian grid. 
Computed in the "box-in-a-star" setu p (for 
details see, e.g., Frevt ag~ et al.l |2012|) . our 
solar model atmosphere was allowed to 
evolve hydrodynamically for several convec- 
tive t urnover times (see iLudwig & Kucinskasl 
I2012L Appendix A, for a definition of different 
time scales). This particular C0 5 BOLD simula- 
tion was also used in the studies of Caffau et 
al. d2008h and B eeck et alJ(l2012h . to which we 
refer the reader for a more detailed description 
of the input microphysics and physical proper- 
ties of the model itself. 

A representative sub-sample of twenty 3D 
model snapshots (i.e. model structures calcu- 
lated at different instants in time) was cho- 
sen out of the relaxed part of the 3D model 
run in order to produce a statistically repre- 
sentative and uncorrected snapshot selection. 
More specifically, snapshots in this sub-sample 
were chosen in such a way that the average 
emerging radiation flux and its standard devi- 
ation would match the corresponding values 
of the entire 3D model run. Similar require- 
ments were applied in the case of the horizon- 
tally averaged vertical mass flux at characteris- 
tic optical depths. The final selection of the 3D 
model snapshots obtained according to these 
criteria was subsequently used in the evalua- 
tion of NLTE effects (Sect. I2.2I ) and, later, in 



the spectral synthesis calculations of the 1 IR 
triplet (Sect. CO). 

The original model snapshots had 140 grid 
points in each horizontal direction and 150 grid 
points on the vertical axis, corresponding to 
a spatial coverage of 5.6 x 5.6 x 2.27 Mm. 
The vertical grid of the model spanned the op- 
tical depth range of -6.7 < log t Ross < 5.5, 
which is sufficient to cover the depths where 
the Oi IR triplet lines form. For the calcula- 
tion of departure coefficients and spectral line 
synthesis, a coarser 3D model atmosphere was 
constructed by choosing every third grid point 
horizontally, thereby reducing horizontal res- 
olution of the model to 47 x 47 grid points. 
We have verified that differences in the spec- 
tral synthesis results obtained with the full and 
reduced model atmospheres, respectively, were 
negligible. 

2.2. NLTE calculations 

Departures from LTE in the line forma- 
tion computations are quantified via three- 
dimensional sets of departure coefficient^] 
which in our stud y were computed w ith the 
NLTE 3D code (seelCavrel et alJl2007t Steffen 
et al"[2012). We have recently made numer- 
ous improvements and generalizations to the 
code in order to adapt it for a wider variety 
of model atoms and astrophysical tasks. More 
specifically, we included IONDIS/OPALAM 
routines, originally part of the Kiel stellar 
atmosphere package. These routines provide 
partition functions for a variety of chemical 
species and are also used for the computa- 
tion of LTE population numbers and contin- 
uous opacities at the wavelengths of the line 
transitions in a given model atom. Since these 
routines are also used in the spectrum synthe- 
sis code Linfor3D (Sect. 12.31 ), this also con- 
tributes towards the self-consistency between 
the NLTE3D and Lin£or3D packages. 

During the first step, the NLTE3D code 
computes photo-ionization and collisional 
rates, where the upward and downward col- 



1 The departure coefficient for atomic level i = 
1 . . . ( max is defined as b t = N iNLTE /N iLTE , where JV,- 
is the level population for the respective case. 
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lisional rates are treated in detailed balance 
and, hence, only the upward rates need to 
be computed explicitly. These rates are kept 
constant throughout all further computations. 
After this step, radiation transfer calculations 
are done for each of the considered line transi- 
tions, £ = 1 . . . /"max, to compute the mean line 
intensity J((x, y, z) at all grid points, where J 
refers to the angle-averaged intensity, and over- 
bar indicates averaging over the line profile. 
Typically, we use 16 inclined plus the verti- 
cal direction for angular averaging, and 37 fre- 
quencies to resolve the line profile. Then, the 
photo-excitation rates are computed using J(, 
and all rates at each grid point of the 3D model 
snapshot are passed to the statistical equilib- 
rium routines that compute NLTE population 
densities according to the prescriptions given 
in Mihalas (1970). Given Jf, a set of linear 
equations needs to be solved for the unknown 
departure coefficients i = 1 . . . / max , which 
is done using standard linear algebra routines, 
independently for each 3D model grid cell. 

Once the new iteration of NLTE population 
numbers is completed, the NLTE line opaci- 
ties are updated and the next iteration is started 
to compute line radiation transfer in order to 
obtain updated photo excitation rates. This cy- 
cle is repeated via the ordinary A-iteration 
scheme until the relative change in the selected 
NLTE line equivalent widths (EWs) becomes 
less than 10~ 3 per iteration (see also Sect. 13. 2\ . 

2.2.1 . The oxygen model atom 

The Grotrian diagram of the oxygen model 
atom used in our study is shown in Fig. Q] The 
model atom consists of 2 max = 16 levels, con- 
nected by 16 bound-free and ^ max = 31 bound- 
bound transitions. Atomic data that describe 
the levels (energies, statistical weights) and ra- 
diative transitions (Einstein coefficients) were 
taken from the NIST database. Photoionization 
cross-sections for the radiative bound-free 
transitions were ado pted from the T opbase of 
the Opacity Project dCunto et al.lll993l and fur- 
ther updates). 

Excitation and ionization due to collisions 
with neutral hydrogen atoms was treated via 
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Fig. 1. Oxygen model atom used in this work. 
Radiatively allowed transitions are marked by 
solid (green) lines, while radiatively forbidden 
transitions that were treated only via electron- 
impact excitation are marked by dashed (blue) 
lines. The forbidden [Oi] line and the Oi IR 
triplet are marked by dash-dotted (red) lines. 
Collisional and radiative bound-free transitions 
were taken into account for each level. 



th e iDrawinl (1 19691) formula in the formulation 
of lLamberidl993h . In the present work, we set 
5h, a parameter that scales the hydrogen col- 
lision rates, to 1.0. Collisional ionization by 
electrons w as treated usin g the classical pre- 
scription of ISeatonl dl962 ). Rate coefficients 
for the excitation by inelastic elect ron colli- 
sions were taken from the work of iBarkleml 
d2007h . allowing us to include the collisional 
transitions between the triplet and quintet sys- 
tems. Finally, we include a resonant charge 
transfer reaction O + H + =± + + H° for 
the gr ound level of oxygen, in th e prescrip- 
tion bv lArnaud & Rothenfluglril985l) . The latter 
process is very efficient in late-type stars and 
therefore ensures that the ground level of O n 
is in LTE with respect to the ground level of 
Oi. 
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For the computation of the photo-ionizing 
radiat ion flux, opacity di stribution functions 
from ICastelli & Kuruczl ((2004) were used, 
while the continuous opacity was calcu- 
lated with the IONDIS/OPALAM routines. 
Background opacities stemming from the 
bound-bound transitions of elements other than 
oxygen are not yet included in the current ver- 
sion of the NLTE3D code. Tests have shown, 
however, that they are not important in the case 
of the oxygen model atom used in this study. 

In the statistical equilibrium calculations, 
we treated the O i ground level (2p 3 P) and the 
upper level of the O i IR triplet (3p 5 P) as sin- 
glets, while, when performing spectrum syn- 
thesis calculations, we have taken the fine- 
splitting into account, assuming that the depar- 
ture coefficients are identical for all sub-levels: 
due to their close energetic values, they should 
be efficiently thermalized and thus experience 
identical sensitivity to NLTE effects. 

2.3. Spectral synthesis calculations 

Spectral synthesis calculations of the Oi 
IR triplet lines were done with Linfor3E0. 
Firstly, Lin£or3D calculates the LTE line 
source function, Slte(v) = B v (Planck's black 
body function) and line opacity, klteM for 
any given line transition. Whenever NLTE 
spectrum synthesis is requested, both the line 
source function and the line opacity are modi- 
fied according to: 

Snlte(v) _ ex P(ff ) ~ 1 

Slte(v) (b low /b up ) exp(^) - 1 ' ' 
and 

^NLTE(V) _ eX P(ff)-(frup/frlow) 

klte(v) ~ blow exp(&)-l ' ( ) 

where b\ ovl and b up are the departure coef- 
ficients of the lower and upper level of the 
transition, respectively. Once the line source 
function and line opacity are computed ei- 
ther in LTE or NLTE, the radiation transfer 
is calculated along 16 inclined plus the ver- 
tical direction, taking into account the differ- 
ential Doppler shifts along each line-of sight. 

2 http://www.aip.de/~mst/Linfor3D/linfor.pdf 



This procedure is repeated for every of the 
twenty 3D model snapshots. Finally, the result- 
ing spectral lines synthesized along each of the 
different directions are combined into a single 
spectral line profile by spatial, temporal, and 
angular averaging. 

3. Results 

3.1. Synthetic Oi IR triplet spectra 

We have conducted NLTE and LTE spectrum 
synthesis calculations of the O i IR triplet lines 
using a 3D solar model atmosphere computed 
with the C0 5 B0LD code. Fig. shows the ob- 
tained synthetic NLTE and LTE Oi IR triplet 
spectrum (identical oxygen abundances were 
used in the LTE and NLTE computations). 
It can be immediately seen that both shape 
and strength of the O i IR triplet lines is sen- 
sitive to NLTE effects: while the line wings 
are very similar in both cases, deeper line 
core is present in NLTE. Hence, in accordance 
with the previous studies, NLTE effects lead to 
higher equivalent widths of the O i IR triplet 
lines and to negative 3D NLTE-LTE abun- 
dance corrections. Such behavior is qualita- 
tively very similar to that reported bv Asplund 
et al. d2004l) or lPereira et al.l(l2009l) . which was 
found to be one of the major factors leading to 
a lower solar oxygen abundance. 

The nature of the NLTE spectral line for- 
mation can be investigated by examining the 
dependence of NLTE effects on the continuum 
intensity across the stellar surface, where high 
and low continuum intensity represents gran- 
ular and intergranular regions, respectively. 
Fig. [3] shows the distribution of the NLTE to 
LTE ratios of local EWs computed from the 
vertical ray (disk-center intensity) and plot- 
ted versus the local relative continuum inten- 
sity. While the plot reveals significant differ- 
ences with respect to the results obtained by 
lAsplund et all (I2004L Fig. 4), who found di- 
minished NLTE EWs in some (tough not all) 
intergranular lanes, our results show a quite 
good quali tative agreeme nt with the computa- 
tions of iKiselmanl (fl993. Fig. 8) (note, how- 
ever, that the adopted Su values were differ- 
ent in all three investigations). Our results indi- 
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Fig. 2. Synthetic 3D NLTE (black line) and 3D LTE (red dots) O i IR triplet lines in the solar spec- 
trum calculated with identical oxygen abundance (disk-integrated flux spectrum). Spectral line 
synthesis calculations were done using twenty 3D model snapshots of a solar model atmosphere 
computed with the C0 5 B0LD code. 



cate that, in NLTE, line strengths are enhanced 
at every position on the disk and NLTE ef- 
fects are more pronounced in the int ergranular 
region s (low continuum intensity). iKiselmanl 
( 1993) suggested that this happens because in 
the intergranular regions the triplet lines form 
higher in the atmosphere where the departures 
from LTE are more severe. 



The origin of the NLTE effects can be 
examined by comparing NLTE and LTE line 
source function s and line opac ities (see, e.g., 
the analysis of Kise lmarJll993h . Fig. [4] shows 
the ratio of NLTE and LTE line source func- 
tions (top panel) and line opacities (bottom 
panel) for the 777.19 nm component of the 
O i IR triplet, as given by Equations ([TJ and 
(0. As was stated previously, the 3p 5 P sub- 
levels share identical departure coefficients, so 
the factor bi ov/ /b up that enters the equations and 
hence the ratios S NLTE /S lte and k nlte /k lte 
are identical for each component of the triplet. 



The plot shows that the source function 
decreases below the local Planck function at 
log tr oss < 0.0. Fig. 3] also contains two spa- 
tially and temporally averaged flux contribu- 
tion function s for the line depression (see 
Magainl ["l986h that were computed at the line 
center (777.1954 nm) and wing (777.15 nm) 
of the line profile. It can be seen that the line 
wings form in a region where Snlte is iden- 
tical or very similar to B v , while the line core 
forms where significant departures from LTE 
may occur. 

On the other hand, the line opacity of the 
O i IR triplet stays close to LTE in the entire 
line forming region. The NLTE/LTE ratio of 
line opacities depends mostly on bi ow , which 
is close to unity in the entire optical depth 
range where the O i IR triplet lines form, while 
the sensitivity to b\ ow /b up ratio is rather small. 
Hence, the influence of NLTE effects on the 
line opacity is weak, and therefore it can be 
inferred that the strengthening of the O i IR 
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Fig. 3. Probability density plot showing the 
NLTE to LTE equivalent width ratio of the 
777.19 nm component of the OI IR triplet, 
plotted against the normalized local continuum 
intensity for all twenty 3D model snapshots 
that were used in the calculations. Both con- 
tinuum intensity and equivalent width refer to 
the disk-center solar spectrum (vertical rays). 
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triplet happens solely due to the diminished 
line source function. According to Fabbian et 
al.(|2009, and references therein), this is caused 
by the photon losses in the line. 

3.2. Numerical issues 

Fig. [5] shows the evolution of the equivalent 
width of the 777.19 nm line (top panel) and 
the relative change of EW (bottom panel) as a 
function of the iteration number. As stated pre- 
viously, the A-iterations are continued until the 
relative change in all NLTE equivalent widths 
of the selected radiative transitions (green lines 
in Fig.[TJ becomes less than 10~ 3 per iteration. 
Even though the triplet line shown in Fig. [5] 
satisfies the convergence criterion already af- 
ter a small number of iterations, certain other 
lines (mainly the resonant UV lines, not shown 
in the figure) may show much slower conver- 
gence and may thus require significantly more 
iterations to obtain the final NLTE population 
densities. 



The top panel of Fig. [5] shows that there is 
a wide range in the number of iterations re- 
quired to arrive at the convergence limit: de- 
pending on the selected 3D model snapshot, it 
may be as low as ~ 30 or as large as ~ 500 
(in the most extreme cases the EWs of O i IR 
triplet lines continue to grow even after ~300 
iterations). One should note, however, that the 
poor convergence properties of the ordinary 
A-iteration scheme pose a well-known prob- 
lem in radiation transfer calculations. To quote 
lHubenvl(l2003h : ' 'it exhibits a pathological be- 
havior in that the solution stabilizes (i.e. rel- 
ative changes of the source function become 
extremely small) long before the correct so- 
lution is reached". All this indicates that a 
more efficient iteration scheme needs to be im- 
plemented into the NLTE3D code in order to 
make the computations more efficient and ac- 
curate. 
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Fig. 5. Top panel: evolution of the disk- 
integrated equivalent width (EW) of the Oi 
777.159 nm line during the iteration process. 
Bottom panel: change per iteration in the EW 
of a given Oi triplet component. The figure 
shows the iteration progress in five individual 
3D model snapshots. Note that in the bottom 
panel data from the individual snapshots merge 
into one single curve. 



4. Conclusions 

We provide a short summary of our first results 
regarding the 3D NLTE spectrum synthesis of 
the O i IR triplet lines, carried out using a 3D 
hydrodynamical solar model atmosphere com- 
puted with the CO 5 BOLD code. Atomic level 
population numbers were calculated using a 
significantly expanded and improved version 
of the NLTE3D code. Our results show that, 
in agreement with previous studies, the Oi 
IR triplet lines are sensitive to NLTE effects, 
which lead to deeper line cores than those ex- 
pected in LTE. This happens mainly because 
in the range of optical depths where the cores 
of the O i IR triplet lines form the line source 
function becomes smaller than the local Planck 
function. A significant deepening of the core, 
along with a net strengthening of the equivalent 



width of the line, leads to negative 3D NLTE- 
LTE abundance corrections. 

The current version of the NLTE3D code 
utilizes an ordinary A-iteration scheme for the 
solution of non-linear statistical equilibrium 
equations. Such approach may be adequate in 
cases where all line transitions are weak and 
thus the radiative rates in the lines are only 
weakly coupled to the level populations (e.g. 
lithium). This, however, is obviously not true 
in the case of oxygen. For a proper treatment of 
this more general situation, the implementation 
of a faster and more reliable iteration scheme 
via some form of accelerated A-iteration (e.g . 
lOlson et"al ] |l986t iRvbicki & Hummel 1 199 11) 
may be needed. This will be amongst the main 
priorities for the further development of the 
NLTE3D code, in order to ensure a higher re- 
liability of the numerical results and a better 
versatility of the code in further applications 
involving more complicated model atoms. 
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